Phase diagram and critical exponents of a Potts gauge glass 
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Abstract 

The two-dimensional g-state Potts model is subjected to a Zq symmetric disorder that allows 
for the existence of a Nishimori line. At </ = 2, this model coincides with the ±J random- 
bond Ising model. For q > 2, apart from the usual pure and zero-temperature fixed points, the 
ferro/paramagnetic phase boundary is controlled by two critical fixed points: a weak disorder point, 
whose universality class is that of the ferromagnetic bond-disordered Potts model, and a strong 
disorder point which generalizes the usual Nishimori point. We numerically study the case q = 3, 
tracing out the phase diagram and precisely determining the critical exponents. The universality 
class of the Nishimori point is inconsistent with percolation on Potts clusters. 

During the last decade, the study of disordered systems has attracted much interest. This is true 
in particular in two dimensions, where the possible types of critical behavior for the corresponding 
pure models can be classified using conformal field theory [Q. Recently, similar classification issues for 
disordered models have been addressed through the study of various random matrix ensembles |^ , but 
many fundamental questions remain open. 

An important category of 2D disordered systems is given by models where the disorder couples to 
the local energy density. Two paradigmatic members of this class are the ±J random-bond Ising model, 
and the g-state ferromagnetic random-bond Potts model. The model to be studied in the present Letter 
can be thought of as an interpolation between these two members; we shall therefore begin by recalling 
some of their basic properties. 

The random-bond Ising model (RBIM) is defined by the energy functional 

"Rising — ^ JijS{Si,Sj) , (1) 

where the sum is over the edges of the square lattice, Si = ±1 are Ising spins, and 6{., .) is the Kroneckcr 
delta function. The random bonds take the values = ±1 according to the probability distribution 

P(Jy) = pd{J,, - 1) + (1 - p)S{J,j + 1) . (2) 

The salient feature of this model is that it marries disorder with frustration, leading to the possibility 
of spin glass order. 

Its phase diagram is generally believed to be as in Fig. [|.a [||. The boundary FP between the 
ferromagnetic and the paramagnetic phases is controlled by three fixed points. The attractive fixed 
points at either end of the phase boundary are respectively the critical point of the pure Ising model 
and a zero-temperature fixed point. Between these two wc find the multicritical point N, intersecting 
the so-called Nishimori line Q 

e^ = {l-p)/p. (3) 

On this line, the replicated version of the model possesses a local Z2 gauge symmetry that, among other 
things, allows for exactly computing the internal energy and for establishing the pairwise equality of 
correlation functions 

[{S.,S..---S,,r^-']^[{S.,S.,---S,,f-], (4) 

where (• • •) denotes the thermal and [• • •] the disorder average. Since the Nishimori line is also invariant 
under Renormalization Group (RG) transformations its intersection N with the FP boundary must 




Figure 1: Phase diagram of the ± J random-bond Ising model (a) and the q> 2 state Potts gauge glass 
(b). 

be a fixed point. However, the widespread belief that the corresponding universality class is of the 
percolation type has recently been refuted on the basis of numerical evidence . 

The other model of special interest to us is the random-bond Potts model (RBPM), which is also 
defined by (0), except that the spins now take q different values, Si = 1,2, q. The most well-studied 
case is that of purely ferromagnetic bonds, such as 

P{J.j) = ls{J^, + Jl) + l^iJ.J + -h). (5) 

with i? = J2/J1 > 1 adjusting the disorder strength. 

In contradistinction to the Nishimori point, the fixed point of this model is situated at weak disorder. 
For q > 2 the disorder is relevant ||], and the corresponding line of fixed points tends to the one of 
the pure Ising model in the limit q ^ 2. As a consequence, the critical exponents can be computed 
perturbatively in a (g — 2)-expansion [Q. According to the RG picture, for g > 2 any small amount of 
disorder should induce a fiow towards the random fixed point. That this is also true for q > 4, where the 
phase transition in the pure model is of the first order, is the content of the Aizenman-Wehr theorem 

§. 

In this Letter we shall consider the model 

n = -Y,5^'i\S,^Sj^,h,), (6) 

where Si = \,2, . . . ,q, and 5^''^ {x) = 1 if x = mod q and zero otherwise. The randomness now takes the 
form of a local "twist" Jij, which is clearly a more severe type of disorder than simple bond randomness. 
The variables Jij are taken from the distribution 

9-1 

P{.h,) = (1 - (<z - l)p)5{J^,) + P E ^^J^i - J) ' (7) 

with < p < l/{q — 1) controlling the strength of the randomness. We shall refer to this model, 
which was originally introduced in Ref. as the Potts Gauge Glass (PGG). The particular form of the 
randomness ensures the existence of a Nishimori line (see below). For q = 2, the PGG reduces to the 
RBIM, and for p = l/g it was studied analytically in It is also connected to the RBPM: To wit, 
when q > 2 the pure Potts model {p = 0) should be unstable to a small amount of randomness, meaning 
that the RG fiow cannot be as indicated on Fig. |l|.a. Instead, we are forced to assume the existence of 
a new fixed point F, intermediary between the pure model and the Nishimori point (see Fig. ^.b). But 
whenever {q — 2), and hence the value of p at F, is sufficiently small, frustration effects are negligible, 
and we should fiow to the same random fixed point as in the RBPM. For reasons of continuity we expect 
this argument to hold true also for higher values of q. 
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The expression of the Nishimori hne was obtained in Ref. but since our notation is shghtly 
different we shah repeat the argument here. We first reexpress the disorder distribution as 

P(Jy) = pe^''*"^'^'^) with K = log {1/p -{q- 1)). (8) 
Consider then the disorder averaged internal energy 




K5'-''\jij) 



-0n 



(9) 



where TV = — l/(g — 1 + e^Y^ , N being the number of sites of the square lattice. Ti is then invariant 
under the gauge transformation Si ^ Si — ai,Jij Jij + Ui — Uj, though P{Jij) is not. Still, E is 
invariant since we sum over all configurations of the disorder. Then, averaging over all the possible gauge 
transformations we get 

Imposing K = (3, there is a remarkable simplification: 

Thus E is regular, and Eq. (||) with K ^ (3 defines the generalized Nishimori line. 
Normalized two-point functions are defined by 

{S,S,) = (g- - 5,)) - l). (12) 

Let us now recall how Eq. can be derived for q ~ 2. We consider n — I for simplicity. Using the 
trivial identities J^') (AS" - Act) = J2l=o '^^'^ - 0'^'''' (Act + and ^^^o - = 1 one readily 

establishes that 

2(5(2) (2(5(2) ( A^) „ 1) (2(5(2) ( A^) - l) , (13) 

and, using the same gauge transformation as before, 

[2(,5(2)(^, - S,)) - 1] = [(2(J(2)(^, - S,)) - in (14) 

This relies crucially on the fact that the above trivial identities generate only two terms, and for general 
q we do not expect simple relations like (||)|^. 

We now turn to our numerical results. Random transfer matrices in the Fortuin-Kasteleyn (FK) 
representation ||ll[ have been a very powerful tool for studying the RBPM ||l^. Unfortunately, the 
random twist variables Jij present in (^) complicate the definition of the FK clusters: only those clusters 
are allowed for which Jy — mod q for any path 7 within the cluster |l^ . It is not obvious how this 
constraint can be generalized to real values of q, and even for integer q keeping track of the necessary 
local information would greatly increase the number of basis states needed. We have therefore found it 
more convenient to write the transfer matrices directly in the spin basis. We work at <; = 3 throughout, 
but expect our conclusions to extend to arbitrary q > 2. 

As we have shown in an earlier publication , the phase diagram can be traced out by investigating 
the effective central charge. To this end we have computed the free energy f'^^ = ^"j^/' on strips of 
various widths L and practically infinite length, M = 10^. The (effective) central charge c can then be 



simple relation for a chiral-type correlator was established in Ref. jgj, but in terms of Eq. Q this does not lead to 
degeneracy in the multiscaling spectrum. 
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p 


L = 3,4 




L = 4, 5 




L = 2,3,4 




L = 3,4, 5 






P 


c 


p 


c 


P 


c 


P 


c 


0.01 


1.0521(5) 


0.76825(3) 


1.0520(5) 


0.78074(9) 


1.0520(5) 


0.79460(6) 


1.0520(5) 


0.7987(3) 


0.02 


1.1061(5) 


0.76874(6) 


1.1061(5) 


0.7815(2) 


1.1061(5) 


0.7953(1) 


1.1061(5) 


0.7998(6) 


0.03 


1.1692(5) 


0.76907(9) 


1.1691(5) 


0.7816(2) 


1.1692(5) 


0.7957(2) 


1.1691(5) 


0.7997(6) 


0.04 


1.244(1) 


0.7685(1) 


1.245(1) 


0.7822(7) 


1.244(1) 


0.7951(3) 


1.245(1) 


0.8020(17) 


0.05 


1.336(1) 


0.7663(3) 


1.337(1) 


0.7799(9) 


1.336(1) 


0.7925(6) 


1.338(1) 


0.7995(25) 


0.06 


1.453(2) 


0.7620(3) 


1.456(2) 


0.7739(11) 


1.454(2) 


0.7882(6) 


1.456(2) 


0.7911(30) 



Table 1: Parametrisation of the ferro/paramagnetic phase boundary. 



obtained as the universal coefficient of the finite-size correction to the free energy for periodic boundary 
conditions |l^ 

+ ^ + (15) 

According to Zamolodchikov's c-theorem [Q, here applied to a non-unitary theory, the effective 
central charge increases along the RG flows and coincides with the (true) central charge at the fixed 
points. The FP boundary (cf. Fig. |l|) can be traced by identifying the maximum of c as a function of T, 
for various fixed values of p. 

Since the randomness is strong, and since the fits to ( |l5|) must be based on at least two different 

sizes L to eliminate the non-universal quantity , we have taken several precautions in order to obtain 

iv) 

small error bars on the /]f . First, for any fixed value of p we use the same realization of the disorder 
for the computations at different values of T. Second, for each strip of length M = 10^ we work in a 
canonical ensemble, meaning that disorder realizations for which the fraction of bonds Jij = J does not 
exactly equal p for each J=l,2,...,g— 1 are discarded. Third, for each strip we average f^jf' over up 
to 10^ independent realizations. 

In Table |l| we show the resulting values of c and the inverse temperature (3 = 1/T at the FP 
boundary. The two-point fits are based directly on (^5|) , while the three-point fits include an additional 
non-universal 1/i* correction ||l^. The existence of an attractive fixed point at p ^ 0.04 with a central 
charge slightly larger than Cpure = 4/5, characterizing the pure 3-state Potts model, is brought out very 
clearly. 

The reader may wonder why data for such small system sizes can possible give any reliable information 
about the thermodynamic limit. Comparison with the pure model {p — 0) shows however that in 
particular the three-point fits converge very rapidly towards the exact result: 03^4 = 0.76803, C4 5 — 
0.78043, C2,3,4 = 0.79431, 03,4,5 = 0.79831 Q. We have extrapolated the data at the fixed point F by 
assuming that for each fit, the relative deviation from the infinite-size result is the same as in the pure 
model. In this way we arrive at the final result 

CF = 0.8025(10), (16) 

which compares favorably with the perturbative result Cpert — |§g§ + 0{q — 2)^ « 0.8026 Q for the 
ferromagnetic RBPM. 

To numerically locate the Nishimori point we measure CcS along the Nishimori line. Since in this case 
p is a function of /? (see Eq. (^ with K = P) we can no longer work in the canonical ensemble of disorder 
realizations. Accordingly our error bars are larger. It is however a big advantage to know the exact 
parametrisation of the Nishimori line, since otherwise we would have had to scan a two-dimensional 
manifold of parameter values [ p8| . 

From the data in Table ^ we conclude that the fixed point N is located at pN — 0.0785(10). Using 
the same extrapolation procedure as above we also estimate 

CN = 0.756(5). (17) 
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p 


L = 3,4 


L = 4,5 


L = 2,3,4 


L = 3,4,5 


0.077 
0.078 
0.079 
0.080 


0.7208(4) 
0.7212(5) 
0.7218(5) 
0.7213(7) 


0.7284(22) 
0.7346(27) 
0.7316(22) 
0.7292(24) 


0.7374(8) 
0.7374(10) 
0.7386(11) 
0.7379(14) 


0.739(6) 
0.754(7) 
0.746(6) 
0.741(6) 



Table 2: Effective central charge along the Nishimori line. 



This is in remarkable agreement with the value of the central charge for the percolation limit in the 
RBPM: c = ^^^^ « 0.7571 |l|. Below we shall return to the question whether the Nishimori point is 
"just" percolation. 

We have also measured magnetic multiscaling exponents r/fe, defined in the plane by 

[{S{x^,y^)S{x2,y2)?''] = {{x2 - x^f + [y^ - yifY""-'^ (18) 

for any integer k. On the semi-infinite cylinder of circumference L, with x G [1, L] and y g] — oo, -l-oo[, 
this reads, using a conformal mapping, 

[{S{x,,y)S{x2.y)r] (x fsin A . (19) 



For a pure system, rjn — nxrj, while for percolation over Potts clusters all rjn coincide. The principal goal 
is here to establish the non-trivial multiscaling at N, rather than to determine the rj^ with extraordinary 
precision. The largest system size employed was L — 12, and we approximate the semi-infinite cylinder 
by taking a length of M — 400L. All runs were averaged over 10^ disorder configurations. 

In Fig. H, we show effective values of r/i{L) along the Nishimori line, for various p close to pN- These 
values were obtained by fitting data for all X2 — a;i = 1, • . • , L/2 to (p^); to judge the systematic error 
due to the inclusion of the smallest Ax = |a;2 — we also display a similar plot for ordinary percolation, 
where ?7porc — 5/24 ~ 0.2083 is known exactly. At the fixed point, r]i(L) must tend to a constant, and we 
conclude that pn — 0.079 — 0.080 with rji = 0.20 — 0.21. Discarding the smallest Ax leads to consistent 
results, but with larger error bars. 

Although our value of rji is consistent with percolation, this scenario can be excluded by considering 
higher moments. E.g. for p — 0.080 and L = 12 we obtain 

?7i = 0.21239(35) 772 = 0.25192(39) 

773 = 0.30824(47) 7/4 = 0.33773(52), (20) 

the corresponding values for p = 0.079 being some 6 % smaller. 

Further evidence against percolation can be obtained by similarly considering the energy-energy 
correlations. In analogy with the RBIM case we associate this with a deviation from N along the 
vertical direction on Fig. |.b. The results for rjl are shown in Fig. ^, and once again we compare with 
the percolation value ?7porc = 2(2 — l/i/perc) = 5/2. In this case, the exponents depend less on the 
precise value of pn, but the finite-size corrections are larger. Extrapolating, we find a value of roughly 
ril = 2.75 — 2.85, rather close to the one obtained for the RBIM Nishimori point = 2.83(2) using a 
similar fit jl^. Discarding data with small Ace leads to larger error bars, but is still consistent with 
ijI ~ 2.85. We have also verified that the energy correlations exhibit genuine multiscaling. 

In conclusion, we have studied a g-statc (Potts-likc) generalization of the ±J random-bond Ising 
model that allows for the definition of a Nishimori line. Apart from a weak disorder fixed point that 
coincides with that of the well-studied random-bond Potts model, the model possesses a strong disorder 
point with multiscaling exponents different from those of percolation. The fixed point structure is remi- 
niscent of that found by S0rensen et al. in the context of a ± J like Potts model, which does however 
not possess the gauge symmetry required for defining a Nishimori line. We believe that it would be 
interesting to study whether the critical points of these two models are indeed identical. Open questions 
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Figure 2: Tyi extracted from (pij| ) with Aa; — 1, . . . , L. We also show the corresponding fit for percolation 
(full line) and the exact value ?7porc — 5/24. 
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Figure 3: extracted from ( p9| ) with Ax = 1, . . . , L. We also show the corresponding fit for percolation 
and the exact value Typ^rc — 5/2- 
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concerning our model include the study of its zero-temperature limit, the possibility of reentrance, and 
of its behavior for q > A. It would also be interesting to examine it using a supersymmetric approach. 

We would like to thank J. Cardy, A. Honecker and P. Pujol for useful discussions. 
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